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ERRATUM 



Due to an error in a Ficrz identity published elsewhere [43] , some of the results presented in this paper are also in 
error. Here we outline which results of this paper are correct and which are incorrect, and then briefly discuss the 
consequences for the latter. Further details can be found in [44]. 

The various Fierz identities we presented in this paper (Eqns. (3,4,7) and all the equations of Appendix A) are all 
correct. However, due to parallel-processing of our efforts, our explicit cross section calculation was performed using 
the Fierz identity given in Okun's textbook. This identity should correctly read 



w = \ A A (F 1A G)r(j A )i, (i) 

A 

t— i ; 

t-H \ in correspondence with our Eq. (A7), but does not. In the published expression of Ref [43], the indices {k,i} are 
incorrectly interchanged with {l,m} on one side of the equation, which is equivalent to exchanging F and G on one 
side of Eq. (1). This transposition is thus not an issue for the usual application of the Fierz identity to 2 — > 2 processes, 
O !■ where F — G. However, it becomes an issue for 2 — > 3 processes. As a consequence, the part of our paper that made 
D ' use of Okun's textbook formula contains errors. The Fierz'ed matrix elements of Eqns. (12-15) are incorrect, as are 
C/5 ! the results that follow from them, namely, Eqns. (16, 22, 24), Eqns. (D1-D3), Eqns. (E6 - E8) and Figures 2,3,4,5. 
• The correct results are that the s-wave contributions from the four matrix elements of Eqns. (12-15) cancel exactly 
in the four Fermi limit where M n M x . Since the four Fermi limit gives a zero s-wave result, there are no interesting 
new figures to be presented in this erratum. However, if the four Fermi limit is not adopted, then the cancellation is 
not complete and an unsuppressed s-wave amplitude results. This means that the interesting parameter space is that 
where M v and M x are comparable. In this regime, one should consider not only the 4 diagrams of Fig. 1, but also 
the two additional diagrams in which the W/Z boson is radiated from the rj propagator. This in turn means that the 
leading s-wave amplitude is one power of (M x /M v ) 2 higher than implied by our errant calculation. In Ref. [44], with 
two additional authors, we present the higher order calculation of the unsuppressed s-wave cross section. 
Many of the important qualitative conclusions of this paper still hold: 

(i) W/Z bremsstrahlung can lift helicity suppressions and thus be the dominant annihilation model, albeit for a 
region of parameter space smaller than we originally proposed. 

(ii) antiprotons produced by the decays of the W and Z gauge bosons can prove lethal for models attempting to 
^ ' produce positrons without overproducing antiprotons. 

iy~j , We thank Ahmad Galea for his signifiant help in identifying the errant formulas. 
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Dark matter annihilation to leptons, XX ~~ > ls necessarily accompanied by electroweak radiative 
corrections, in which a W or Z boson is radiated from a final state particle. Given that the W and 
Z gauge bosons decay dominantly via hadronic channels, it is thus impossible to produce final 
state leptons without accompanying protons, antiprotons, and gamma rays. Significantly, while 
many dark matter models feature a helicity suppressed annihilation rate to fermions, radiating a 
massive gauge boson from a final state fermion removes this helicity suppression, such that the 
branching ratios Br(£vW), Br(£ + g~ Z), and H>x(vvZ) dominate over Br(££). W/Z- bremsstrahlung 
thus allows indirect detection of many WIMP models that would otherwise be helicity-suppressed, or 
v 2 suppressed. Antiprotons and even anti-deuterons become consequential final state particles. This 
is an important result for future DM searches. We discuss the implications of W/Z-bremsstrahlung 
for "leptonic" DM models which aim to fit recent cosmic ray positron and antiproton data. 



PACS numbers: 95.35.+d, 12.15.Lk, 95.85.Ry 

I. INTRODUCTION 

An abundance of cosmological and astrophysical evidence 
attests to the existence of dark matter (DM) , whose pres- 
ence is inferred via its gravitational influence [1-3]. How- 
ever, the fundamental particle properties of DM remain 
essentially unknown. One important means of probing 
DM's particle nature is via indirect detection, whereby 
we search for products of DM annihilation (or decay) em- 
anating from regions of DM concentration in the Universe 
today. 

The dark matter annihilation cross section is often 
parametrized as (v<ja) = a + bv 2 + • • • , where (v<ja) is 
the thermally-averaged annihilation cross section. The 
constant a comes from s-wave annihilation, while the ve- 
locity suppressed bv 2 term receives both s-wave and p- 
wave contributions; the L th partial wave contribution to 
the annihilation rate is suppressed as v 2L . Given that 
v <~ lCT 3 c in galactic halos, even the p-wave contribution 
is highly suppressed and thus only the s-wave contribu- 
tion is expected to be significant in the Universe today. 
However, in many DM models the s-wave annihilation 
into a fermion pair XX ^ ff 1S helicity suppressed by a 
factor (m,f/M x ) 2 (only — > it modes remain of interest, 
and then only for a certain range of x mass). 

When computing DM annihilation signals, it is nor- 
mally assumed that only the lowest order tree-level pro- 
cesses make a significant contribution. However, there 
are important exceptions to this statement. Dark matter 
annihilation into charged particles, XX If, 1S neces- 
sarily accompanied by the internal bremsstrahlung pro- 
cess XX f fli where the photon may be radiated from 
one of the external particle legs (final state radiation, 
FSR) or, possibly, from a virtual propagator (virtual 
internal bremsstrahlung, VIB). On the face of it, the 
radiative rate is down by the usual QED coupling fac- 



tor of a J ~ 500. However, and significantly, photon 
bremsstrahlung can lift the helicity suppression of the s- 
wave process [4], which more than compensates for the 
extra coupling factor. Such a striking enhancement can 
arise when a symmetry of the initial state \X is satis- 
fied by the three body final state ffj, but not by the 
two body final state //. For bremsstrahlung of photons, 
only VIB is effective in lifting the helicity suppression, 
as FSR is dominated by soft or collinear photons (such 
that the two and three body final states have the same 
symmetry properties) as discussed in Ref. [5]. 

In this paper we examine electroweak 
bremsstrahlung [6-12], i.e., bremsstrahlung of Z or 
W ± electroweak gauge bosons to produce ffZ and 
tvW final states. The virtue for W/Z bremsstrahlung 
to lift initial-state velocity and final-state helicity sup- 
pressions, alluded to in [8, 11], has not been previously 
explored. We show that VT/Z-bremsstrahlung can also 
lift suppression and become the dominant annihilation 
channel. Thus, W/Z bremsstrahlung allows indirect 
detection of many WIMP models that would otherwise 
be helicity-suppressed, or v 2 suppressed. This is an 
important result for future DM searches. 

There are a number of important distinctions be- 
tween electromagnetic (EM) and electroweak (EW) 
bremsstrahlung. An obvious one is that EM 
bremsstrahlung produces just photons, whereas EW 
bremsstrahlung and subsequent decay of the gauge 
bosons leads to leptons, hadrons and gamma rays, of- 
fering correlated "multi-messenger" signals for indirect 
dark matter searches. Another distinction is that W/Z- 
bremsstrahlung from final state particles (FSR) is suf- 
ficient to lift a suppression. This is due to the nonzero 
gauge boson masses, and the coupling of the gauge bosons 
to the non-conserved axial current which leads to a differ- 
ent form for the polarization sum than in the case of the 
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photon (or a gluon in the similar QCD process). In con- 
trast, for the EM process, VIB is required for the photon 
to lift a suppression. Because an additional propagator 
appears for VIB, suppression-lifting EM bremsstrahlung 
is itself suppressed by an additional factor of M 2 /M 2 rel- 
ative to electroweak's FSR, where M v is the mass of the 
internal exchange-particle. Only in the event of a near- 
degeneracy M x ~ M v is this relative suppression of EM 
bremsstrahlung negligible. 

DM annihilation to charged leptons has been the sub- 
ject of much recent attention, due to recently measured 
cosmic ray anomalies which point to an excess of cos- 
mic ray positrons above those that may be attributed to 
conventional astrophysical processes. PAMELA has ob- 
served a sharp excess in the e + / (e~ +e + ) fraction at ener- 
gies beyond approximately 10 GeV [13], without a corre- 
sponding excess in the antiproton/proton data [14, 15], 
while Fermi and HESS have reported more modest ex- 
cesses in the (e~ +e + ) flux at energies of order 1 TeV [16]. 
These signals have led to a re-examination of positron 
production in nearby pulsars [17], emission from super- 
nova remnants [18], acceleration of e + e~ in cosmic ray 
sources [19], and propagation in conventional cosmic ray 
models [20]. As an alternative to these astrophysical 
mechanisms, it has also been proposed that the excess 
e + and e~ are produced via dark matter annihilation 
in the Galactic halo, with an abundance of DM models 
proposed to accomplish this end. A recent overview of 
e ± -excess data and possible interpretations is available 
in [21]. 

However, some of the most popular models suffer from 
helicity or w 2 -suppression. A prototypical example of 
suppressed production of Standard Model (SM) fermion 
pairs is provided by supersymmetry: Majorana neutrali- 
nos annihilate into a pair of SM fermions via t- and it- 
channel exchange of S , f/(2)-doublet sfcrmions. To over- 
come the suppression, proponents of these models have 
invoked large "boost" factors. These boost factors may 
be astrophysical in origin, as with postulated local over- 
densities of dark matter, or they may arise from particle 
physics, as with the Sommerfeld enhancement that arises 
from light scalar exchange between dark matter particles. 
Although not ruled out, these factors do seem to be a con- 
trivance designed to overcome the innate suppression. 

A further problem with suppressed models is the 
overproduction of antiprotons from unsuppressed W/Z 
bremsstrahlung. Given that hadronic decay modes of 
the W and Z bosons will lead to significant numbers of 
both antiprotons and gamma rays, this will impact the 
viability of models that might otherwise have explained 
the observed positron excess. Even in models which do 
not feature a suppression, the W/Z-bremsstrahlung has 
important phenomenological consequences, as the decay 
products of the gauge bosons make a pure leptonic e + e~ 
signal impossible [11]. 

In Section II we discuss the circumstances under which 
dark matter annihilation may be suppressed, and in Sec- 
tion III explain how W/Z bremsstrahlung is able to cir- 



cumvent such a suppression. In Section IV we consider 
a representative model, and explicitly calculate the cross 
sections for both the lowest order annihilation process, 
and for the W/Z bremsstrahlung process. We discuss 
implications of these results in Section V. Calculational 
details are collected in five Appendices. 



II. UNDERSTANDING SUPPRESSION USING 
FIERZ TRANSFORMATIONS 

In this section we describe the origin of v 2 and helicity 
suppressions. We shall make use of Fierz transformation 
and partial wave decomposition to determine under what 
circumstances these suppressions will or will not arise. 

Dark matter candidates may be scalar, fermionic, or 
vector in nature; if fermionic, they may be either Dirac 
or Majorana. Permissible annihilation models include s-, 
t-, and u-channel exchanges of a new particle, and the 
various possibilities are listed in Refs. [22-24]. In every 
case, it is useful to classify the partial waves available to 
the decay process, and to analyze the dependence on the 
mass of the SM particle-pair in the final state. In this 
article, we focus on fermionic Majorana dark matter. 

For fermionic dark matter, the natural projection of 
2^2 processes into partial waves makes use of the 
Fierz transformation. In the next subsection we con- 
sider DM annihilation via the process XX ~^ If: an d 
explain the use of Fierz transforms to convert the matrix 
elements for i/w-channel annihilation, which are of the 
form (x TaI){1 Tbx)j_^° a sum °f s-channel amplitudes 
of the form (% Tix)(l T2I). In the following subsection 
we then categorize the Fierzed s-channel amplitudes into 
partial waves and fcrmion-pair spin states, which deter- 
mines whether the amplitudes are velocity suppressed, 
mass-suppressed, or unsuppressed. In the third and final 
subsection, we put our findings together to determine 
which class of models will have a suppressed 2 — > 2 anni- 
hilation. We show that in a certain popular class of sup- 
pressed models, the 2^3 TV/Z-bremsstrahlung process 
is unsuppressed, and in fact dominant for 2 M x > My/. 
We will find in Section III that a generalization of the 
Fierz transformation offers useful insight into the non- 
suppression of the 2 — > 3 process. 



A. Fierz Transformations in the Chiral Basis 

Helicity projection operators are essential in chiral 
gauge theories, so it is worth considering the reformu- 
lation of Fierz transformations in the chiral basis [25]. 
(A discussion of standard Fierz transformations may be 
found in, e.g. Ref. [26].) We place hats above the gener- 
alized Dirac matrices constituting the chiral basis. These 
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matrices are 

{f B } = {P R , P L , P R -f, P L r, , and 

{Tb} = {Pr, Pl, Pl-Yv, Pr^, \ <W . ( 2 ) 

where Pr = |(1 + 75) and Pl = |(1 — 75) are the usual 
helicity projectors. Notice that the dual of Prj^ is PlJ^, 
and the dual of PlY is PrJ^- The tensor matrices in this 
basis contain factors of \: T T = ■^o lxv and IY = 
These facts result from the orthogonality and normaliza- 
tion properties of the chiral basis and its dual, as ex- 



plained in detail in Appendix A. 

Using completeness of the basis (see Appendix A) , one 
arrives at a master formula which expands the outer 
product of two chiral matrices in terms of their Fierzed 
forms: 

(f D ) [f E ] = l -Tr [t D f c f E t B ] (f B ] [f c ) , (3) 

where the parentheses symbols are a convenient short- 
hand for matrix indices [27] (see the appendix for de- 
tails). Evaluating the trace in Eq. (3) leads to the Fierz 
transformation matrix in the chiral-basis: 



/ (Pr) [Pr] \ 
(Pl) [Pl] 

(f ) [f ] 
(75 f) [f] 
(Pr) [Pl] 
(P R r ) [PliA 
(Pl) [Pr] 

(p L -n [Pria 
(p R r) [Pr^\ 
V (p L r) [PliA J 



(2 1 1 
2 1-1 
6 6-20 
6-602 
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/ (Pr) [Pr) \ 
(Pl] [Pl) 

(f ] [f ) 
(75 T] [T) 
(Pr] [Pl) 

(p R r] \plih) 

(Pl] [Pr) 
(P L r] [P«7m) 
(PrY] [Pr1») 
V (PlY\ [Pl7m) / 



Non-explicit matrix elements in (4) are zero, and we have 

introduced a shorthand T for either T T = ^cr MI/ or Ft — 
1 

2 C 7" y- 

The importance of this transformation for us is that 
it converts ^-channel and w-channel exchange graphs into 
s-channel form, for which it is straightforward to eval- 
uate the partial waves. The block-diagonal structures, 
delineated with horizontal and vertical lines, show that 
"mixing" occurs only within the subsets {Pr<S>Pr, Pl <8> 
P L , t®f, -y b f®f}, and {P R ®P L , Pr1»®Pl1»}- The 
Fierz transform matrix is idempotent, meaning its square 
is equal to the identity matrix. This follows from the 
fact that two Fierz rearrangements return the process to 
its initial ordering. A consequence of the block-diagonal 
form is that each sub-block is itself idempotent. 

In Eq. (4) we have included one non-member of the ba- 
sis set, namely 75 T; it is connected to T via the relation 

75 ® = -c 1 * p a afi . (5) 

Explicit use of 75 T in Eq. (4) is an efficient way to express 
the chiral Fierz transformation. 

So far we have not used the qualifier in the assumption, 
that the dark matter is Majorana. Majorana particles 
are invariants under charge conjugation C, which implies 
that vector and tensor bilinears are disallowed. Another 
way of understanding this is to note that interchanging 
the two identical Majorana particles in a t-channel dia- 



gram generates an accompanying w-channel diagram with 
a relative minus sign (from fermion anticommutation). 
When Fierzed, these two amplitudes cancel for V and 
T couplings (exactly so in the Four-Fermi limit where 
the differing momenta in the t- and u-channel propaga- 
tors can be ignored - refer to Appendix B for details). 
We must thus drop V and T couplings appearing in the 
Fierzed bilinears of the x-current. 



B. Origin of v 2 and Helicity Suppressions 

One can use partial wave decomposition (see. e.g., 
the textbooks [28-30] , or the convenient summary in the 
Appendix of [8]) to expand the scattering amplitudes as 
a sum of angular momentum components. Partial waves 
do not interfere, and the L th partial wave contribution 
to the total cross section av is proportional to v 2L . The 
annihilating \ particles are very non-relativistic today, so 
an unsuppressed s-wave (L = 0), if present, will dominate 
the annihilation cross section. The DM virial velocity 
within our Galaxy is about 10~ 3 (in units of c), leading 
to a suppression of v 2 <~ 10~ 6 for p-wave processes. 

On the other hand, the SM fermions produced in the 
2^2 annihilation are highly relativistic (except possi- 
bly for tt production). For many annihilation channels 
the spin state of the fermion pair gives rise to a helic- 
ity suppression by a factor of (m;/M x ) 2 , where mi is the 
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fermion mass. 

Unfortunately, many popular models for annihilation 
of Majorana dark matter to charged leptons are sub- 
ject to one or more of these two suppressions, the v 
and/or (mi/M x ) 2 suppressions. This includes some of 
the models proposed to accommodate the positron and 
e + e~ excesses observed in PAMELA, Fermi-LAT, and 
HESS data. In Section III, we show that in the class of 
models which have suppressed rates for XX ~^ , the 
2^3 graph obtained by adding a radiative W ± or Z to 
the final state particles of the 2^2 graph becomes dom- 
inant. The radiated W's and Z's will decay to, among 
other particles, antiprotons. Since an excess generation 
of antiprotons is not observed by PAMELA, this class of 
models is ruled out by the present work. 

Consider products of s-channel bilinears of the form 
(x Tix)G r 2 Z)- To further address the question of which 
products of currents are suppressed and which are not, we 
may set v 2 to zero in the x-current, and m 2 to zero in the 
lcpton current, and ask whether the product of currents 
is suppressed. If the product of currents is non-zero in 
this limit, the corresponding amplitude is unsuppressed. 
In Table I we give the results for the product of all stan- 
dard Dirac bilinears. (The derivation of these results is 
outlined in Appendix C.) Suppressed bilinears enter this 
table as zeroes. 1 

One can read across rows of this table to discover that 
the only unsuppressed s-channel products of bilinears for 
the 2 — > 2 process are those of the pseudo-scalar, vector, 
and tensor. (For completeness, we also show results for 
the pseudo-tensor bilinears, although the pseudo-tensor 
is not independent of the tensor, as a result of Eq. 5.) For 
Majorana dark matter, the vector and tensor bilinears 
are disallowed by charge-conjugation arguments and one 
is left with just the unsuppressed pseudo-scalar. 

C. Class of Models for which \X ~^ Annihilation 
is Suppressed 

We now put the results of the previous two subsec- 
tions together to explain which class of models have a 
v 2 and/or (me/M x ) 2 suppressed 2 — > 2 annihilation. We 
have seen that, for Majorana DM, s-channel annihilation 
with a P coupling is unsuppressed, while S and A con- 
tributions are suppressed (and V and T forbidden). Let 
us now consider i-channel or w-channel processes. 

Any i-channel or u-channel diagram that Fierz's to an 
s-channel form containing a pseudoscalar coupling will 



It is seen that the only bilinears in the table without velocity- 
suppression are those of the pseudo-scalar, the three-vector part 
of the vector, the zero th component of the axial vector, and the 
time-space part of the tensor (or equivalently, the space-space 
part of the pseudotensor). It is also seen that the only bilinears 
without fermion mass-suppression are the scalar, pseudoscalar, 
three- vector parts of the vector and axial vector, and the tensor. 



have an unsuppressed L — s-wave amplitude. From 
the matrix in Eq. (4), one deduces that such will be the 
case for any t- or w-channel current product on the left 
side which finds a contribution in the 1 st , 2 nd , 5 th , or 
7 th columns of the right side. This constitutes the t- 
or u-channel tensor, same-chirality scalar, and opposite 
chirality vector products (rows 1 through 4, and 6 and 
8 on the left). On the other hand, the t- or u-channel 
opposite chirality scalars or same-chirality vectors (rows 
5, 7, 9, and 10 on the left) do not contain a pseudoscalar 
coupling after Fierzing to s-channel form. Rather, it is 
the suppressed axial-vector and vector (Dirac fermions 
only) that appears. 

Interestingly, a class of the most popular models for 
fermionic dark matter annihilation to charged leptons, 
fall into this latter, suppressed, category. It is precisely 
the opposite-chirality t- or u-channel scalar exchange that 
appears in these models, an explicit example of which will 
be discussed below. Thus it is rows 5 and 7 in Eq. (4) 
that categorize the model we will analyze. After Fierzing 
to s-channel form, it is seen that the Dirac bilinears are 
opposite-chirality vectors (i.e., ^ or i). Dropping the 
vector term from the x-current we see that the 2 — > 2 
process couples an axial vector x-current to a relativis- 
tic SM fermion-current which is an equal mixture of A 
and V. Accordingly, this model has an s-wave amplitude 
occurring only in the L = 0, J — I, S— 1 channel, with 
the spin flip from S = to S = 1 (or equivalently, the 
mismatch between zero net chirality and one unit of he- 
licity) costing a fermion mass-insertion and a (m//M x ) 2 
suppression in the rate. 

Let us pause to explain why this t- or u-channel scalar 
exchange with opposite fermion chiralitics at the vertices 
is so common. It follows from a single popular assump- 
tion, namely that the dark matter is a gauge-singlet Ma- 
jorana fermion. As a consequence of this assumption, 
annihilation to SM fermions, which are SU(2) doublets 
or singlets, requires either an s-channel singlet boson or a 
t- or it-channel singlet or doublet scalar that couples to \- 
f. In the first instance, there is no symmetry to forbid a 
new force between SM fermions, a disfavored possibility. 
In the second instance, unitarity fixes the second vertex 
as the hermitian adjoint of the first. Since the fermions of 
the SM are left-chiral doublets and right-chiral singlets, 
one gets chiral-opposites for the two vertices of the t- or 
M-channel. 

Supersymmetry provides an analog of such a model. In 
this case the dark matter consists of Majorana neutrali- 
nos, which annihilate to SM fermions via the exchange 
of ("right"- and "left" -handed) SC/(2)-doublet slepton 
fields. In fact, the implementation in 1983 of supersym- 
metric photinos as dark matter provided the first explicit 
calculation of s-wave suppressed Majorana dark mat- 
ter [31]. However, the class of models described above 
is more general than the class of supersymmetric models. 

To illustrate our arguments, we choose a simple ex- 
ample of the class of model under discussion. This is 
provided by the leptophilic model proposed in Ref. [32] 
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TABLE I. Extreme non-relativistic and extreme relativistic limits for s-channel bilincars. In order for a term with an initial- 
state DM bilinear and a final-state lepton bilinear to remain unsuppressed, the DM bilinear must have a non-zero entry in the 
appropriate cell of the "v = limit" columns, and the lepton bilinear must have a non-zero term in the appropriate cell of the 
"M = limit" columns. Otherwise, the term is suppressed. (The tensor and pseudo-tensor are not independent, but rather 
are related by 75 — i e f" /a ' 3 cr a/3 .) We recall that antiparallel spinors correspond to parallel particle spins (and antiparallel 
particle helicities for the M = current), and vice versa. Amplitudes are shown {ot uF D v = \vTr>u]*. The two- fold ± 
ambiguities reflect the two-fold spin assignments for parallel spins, and separately for antiparallel spins. 



by Cao, Ma and Shaughnessy. Here the DM consists of 
a gauge-singlct Majorana fermion \ which annihilates to 
leptons via the SU (2)-invariant interaction term 

/ (ut-) L e ^ j£ j X + h.c. = f{v L rf - ^ry+) X + h.c. 

(6) 

where / is a coupling constant, e is the 2x2 antisym- 
metric matrix, and (rj + , rf) form the new SU(2) doublet 
scalar which mediates the annihilation. (This model was 
originally discussed in Ref. [33] , and an expanded discus- 
sion of its cosmology may be found in Ref. [34].) 

As discussed above, the u- and ^-channel ampli- 
tudes for DM annihilation to leptons, of the form 
(X-PtO (IPrx), become pure (xPlI^x) (IPrJ^I) under 
the chiral Fierz transformation. The product of the Ma- 
jorana and fermion bilinears then leads to an AA term 
and an AV term. However, reference to Table I shows 
that neither of these terms leads to an unsuppressed am- 
plitude: in all cases, either the lepton bilinear is sup- 
pressed by mg, the DM bilinear by u, or both are sup- 
pressed. Thus, Majorana DM annihilation to a lepton 
pair is suppressed in this model, in accordance with the 
explicit calculation in Ref. [32]. 



tor ~ (M x /Mvi/) 2 /87r 2 relative to 2-body massless phase 
space, which exceeds unity for dark matter masses ex- 
ceeding ~TeV. 2 More importantly, the v 2 suppression 
for Majorana annihilation to 2-body final states will be 
lifted by the 3-body VF-bremsstrahlung process. In Sec- 
tion IV we show, by explicit calculation, that the 2^3 
radiative process that leads to antiprotons dominates for 
any M x that allows the W to be produced on-shell, i.e., 
for 2M X > M w . 

The next inevitable question is "Why is the radiative 
2^3 process unsuppressed?" To answer this question, 
we invoke a more general Fierz rearrangement applicable 
to 2 — > N processes, N > Z. The relevant equation, 
derived in Appendix A states that any 4x4 matrices X 
and Y may be expressed as 

(X) [Y]= (XI) [1Y] 

= J(xr s Y] [r B ) 

= lTr[Xr B Yr c ] [T c ] [Y B ), (7) 

where the Dirac matrices here are taken in the standard 
basis defined in Eq. (Al). 

From Table (I) we see that setting T c to 75 7 , the 
only structure available to a non-relativistic Majorana 



III. LIFTING THE SUPPRESSION 



Allowing the lepton bilinear to radiate a W or Z bo- 
son (as shown in Fig. (1)) does yield an unsuppressed 
amplitude. In the rate, there will be the usual radia- 
tive suppression factor of ^ ~ 10~ 3 . But, this will 
be partially compensated by a 3-body phase space fac- 



2 When M 2 2> M^y , the rate for single W production is dominated 
by infrared and collincar divergences, leading to a suppressed fac- 

tOT (3%) " 3111 (3%) I 10 ' 35 1 inStCad ° f ° Ur 
Moreover, the rate for multiple production of VK's becomes so 
large that resummation techniques are necessary. 
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x(fci) 



(A) 



e (gi) v(pi) 
— > 1 >- 



X(ki) 



' W~ (Q) 



(B) 




e (<?i) ^(pi) 
— » 1 >- 



»W- (Q) 



x(fc) 



f(pl) 
— > — 



x(fci) 



(Q) 




" (92) e+ (p 2 ) X (fa) 



— > — 



' W~ (Q) 



v (q 2 ) e+ (p 2 ) 



(D) 



FIG. 1. i-channel (A and C) and u-channel (B and D) contributions to \X e + vW~ . Emission from the scalar propagator 
is not included, as it is suppressed by 1/M V . Note that all fermion momenta flow with the arrow except so qi = pi + Q, 
qi = —P2 — Q. 



current other than the pseudoscalar, and Tb to either 7^ 
or 757 J , provides an unsuppressed product of the Majo- 
rana dark matter bilinear and the charged lepton bilinear. 
Moreover, for the W/Z-bremsstrahlung process, X and 
Y in the general Fierz equation are the un-Fierzed cou- 
plings Pl and q~ 2 Pr^Pl ft respectively. So we will have 
shown that the radiative process is unsuppressed if we 
can show that q~ 2 Tr [Pl (7 j or 757^) Pr$Pl $ 75 70] is 
unsuppressed. This trace reduces to q~ 2 Tr [Pr 70 7 J ft. 
The expansion of this trace as scalar products contains 
terms such as q ■ e J and (e x q)i , which are nonzero and 
unsuppressed by fermion masses. Thus, the 2 — > 3 pro- 
cess contains an unsuppressed s-wave amplitude. 

Physically, the un-suppression works because the 
gauge boson carries away a unit of angular momentum, 
allowing a fermion spin-flip such that there is no longer a 
mismatch between the chirality of the leptons and their 
allowed two-particle spin state. 

One may ask why emission of a gamma-ray rather 
than a W/Z boson is less effectual. It has been known 
for some time [4, 5] that gamma-ray emission in the fi- 
nal state does produce an unsuppressed s-wave contribu- 
tion, but at second order rather than lowest order in the 
inverse mass-squared M~ 2 of the t- and u-channel ex- 
change particle(s). The reason is that gamma-ray emis- 
sion from the final state fermions (FSR) is dominated 
by infra-red and collinear singularities, each of which 
puts the intermediate lepton on-shell (virtuality q 2 — > 0) . 
Including the q~ 4 from the squared propagator in the 
phase space integral (see Eq. (21)), one gets the factor 



Jm 2 lH s _ I 2 ) (<? 2 ~ M v)i where M v is the mass of the 
radiated boson (photon or W or Z). For a gamma-ray, 
with My = 0, one readily sees the infra-red and collinear 

singularities in J Q . An on-shell particle is observable, 
so the spin states of the q 2 — > intermediate fermion do 
not interfere. Thus, as q 2 — > 0, the trace for gamma 
emission, Tr [70 j j 0] = Tr [70 j j (P| + Pl)$ ft goes 
over to Tr [ 7o 7 3 Pr] Tr [P R f ft + Tr [ 7o ^ P L ] Tr [P L $ ft . 
The first trace in each term of this sum vanishes. Con- 
sequently, the gamma-emission amplitude remains sup- 
pressed at order M~ 2 . However, at order M~ 4 , the 
gamma-ray may be emitted from the internal particle 
rj (VIB). For VIB, phase space does not favor q 2 = 0, 
and an unsuppressed amplitude results. 

The emission of a massive W (or Z) boson contrasts 
significantly from the emission of a massless photon. 
With the W emission, the relevant phase space integral 
over virtuality q 2 is ^f-(s-q 2 ) {q 2 -M 2 v ). The min- 
imum virtuality of the intermediate fermion is q 2 — , 
and the mean virtuality for s 3> is greater again by 
the factor 2 ln(s/M^). With no infra-red or collinear sin- 
gularities for WyZ-emission, an unsuppressed amplitude 
results already at order M~ 2 . 

Before looking at an explicit example in which elec- 
troweak bremsstrahlung is seen to lift a suppression, we 
pause to summarize some important facts for the 2^2 
annihilation process: 

• Fierz transformation is used to re-express t- and u- 
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channel amplitudes of the form (x t A l)(l f bx) as a 
sum of s-channel (not to be confused with s-wave) 
amplitudes of the form (x T c x)(l f^Z). 

• For Majorana dark matter, only S, P, and A s- 
channel bilinears are allowed, with the V and T 
bilincars forbidden by the self-conjugate properties 
of Majorana particles. 

• Considering the product of an s-channel %-current 
with an s-channel fermion-current, we find that the 
pseudo-scalar is the only member of the set (S, 
P, A) which is unsuppressed. The other combi- 
nations are either helicity (mg/m x ) or velocity (v) 
suppressed. 

• The annihilation process XX ~^ ^ Yia t~ an d u- 
channel exchange of a scalar is suppressed. Im- 
portantly, electroweak bremsstrahlung lifts this 
suppression at lowest order in the propagator 
mass-squared (M~ 2 in amplitude), whereas photon 
bremsstrahlung lifts the suppression at the next or- 
der (M~ 4 in amplitude). 

Amplification of the latter remark is the purpose of 
this paper. 



interference term. Alternatively, one may Fierz trans- 
form the fermion bilinears in the two contributing am- 
plitudes. The relative minus sign is compensated by the 
special Majorana minus sign described in Eq. (B2). Ref- 
erence to Eq. (4) then shows that one gets {Pl) [Pr] — > 
5 {PlI^} [Pr1h) x 2, where the final factor of 2 counts the 
two contributing amplitudes, which are identical in the 
four-fermi limit M 2 » t and u. We are left with just one 

amplitude, j^\v{k2){\^)v{p2)] [«(pi)Pl7a» u (P2)]. The 
surviving Dirac structure for the Majorana current is 
pure axial vector, since the vector (and tensor) part of a 
Majorana current vanishes. With just a single product 
of bilinears, the remaining part of the 2—^2 calculation 
is straightforward. One arrives at 

in agreement with the four-fermi, mg = limit of Eq. (8). 
Here, the helicity suppression of the s-wave amplitude, 
proportional to a helicity flip, in turn proportional to a 
mass insertion, is manifest. 



B. W Emission and Unsuppressed S- wave 



IV. EXPLICIT CALCULATION OF 
SUPPRESSION-LIFTING WITH 
ELECTROWEAK BREMSSTRAHLUNG 

To explicitly demonstrate that emission of a or 
Z boson does lift helicity suppression, we calculate the 

cross section for XX ~^ ° T vW ± below in the leptophilic 
model of Rcf [32] . The interaction term for this model is 
that given above in Eq. (6). 



A. Example of Helicity-Suppressed Rate 

In the model of ref . [32] , the cross section for the 2^2 
process XX c + e~ or vv with Majorana DM is given as 



v a 



f 4 v 2 r 2 

24ttM 2 



(1 -2r + 2r 2 ), 



(8) 



where mi ~ and M„± = M„o have been assumed, and 
r = M 2 /(M 2 + M 2 ). The suppressions discussed in Sec- 
tion II are apparent in Eq. (8). The helicity suppressed 
s-wave term is absent in the mi = limit, and thus only 
the w 2 -suppressed term remains. 

This 2^2 cross-section can be calculated by inclusion 
of two Feynman diagrams, a i-channcl exchange of 77 and 
the associated ii-channel exchange obtained by crossing 
the Majorana particles. The relative sign between the 
graphs is negative, due to the fermion exchange. Sum- 
ming and squaring, one has three terms including the 



We now turn to the calculation of the cross section 
for the process XX ~^ e + vW~ (equal to that for XX ~^ 
e~i>W + ). The four contributing Feynman diagrams are 
shown in Fig. 1. Note that we consider bremsstrahlung 
only from the final state particles (FSR), and neglect 
emission from the virtual scalar (VIB). Strictly speak- 
ing, the distinction between FSR and VIB is somewhat 
artificial in the sense that the partition depends upon 
the choice of gauge. However, we shall work in unitary 
gauge, in which emission from the internal line is sup- 
pressed by a further power of M 2 due to the additional 
scalar propagator; consequently, we expect our results to 
be valid to order M~ 2 in amplitude, i.e. order M~ 4 in 
rate. 

We retain the assumptions m; ~ and M v ± — M^a. 
The matrix element for the top-left diagram is 



M A = 



Wf 2 



1 



-(v{k 2 )PLv{ P 2)) 



V2q 2 ti - M 2 
xfuip^PL^uihfje®. 



(10) 



where i-^i^Ph is the coupling at the IvW vertex, and 
t\,t2,u\,U2 are the standard Mandelstam variables, 



h = (fci - qi) 2 = (P2 - k 2 f 
*2 = (fci - Pi? = (-92 - k 2 f 
Ml = (k 2 - qi) 2 = (p 2 - fcl) 2 
U 2 = (k 2 -pi) 2 = {-q 2 - fcl) 2 - 



(11) 



Upon applying Eq (7) to Fierz transform the matrix el- 
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ement, we obtain 

wf 2 1 



V2qf h -M% ^ 4 



+ 



^gf 2 



^2q\ h -M* "4 
x (Hhh 5 jMki))(u(Pi)PLl a rfiv( P2 )^ .(12) 

The first two terms after the first equality are zero due 
to the helicity projection operators, leaving only an axial 
vector term. (Vector and tensor x _ bilinears have been 
omitted, as they will cancel between u and t channel di- 
agrams in the heavy M v limit, as discussed above.) Note 
that although this matrix element resembles that of an s- 
channcl annihilation process, the 7 matrices in the lepton 



bilinear would be in a different order for a true s-channcl 
annihilation process involving VK/Z-bremsstrahlung from 
one of the final state leptons. 

Similarly, the matrix element for the top-right diagram 
can be written as 



Mb = 



-*g.f 2 1 1 



x(u( Pl )P Ll a r0iv(P2))e^ (13) 



and those for the bottom diagrams, 



Mi 



M 



V2 q lt 2 -Mf4^ k2)w{h) ) 
x(«(pi)PL^7 M 7 a ^b2))^, (14) 



igf 2 1 1 



— ^2-(v(fc 2 ) 75 7aU(/ci)) 



V2q% u 2 ^ 

x(u(Pi)p L ^rrv(P2))tf- (15) 

Performing the sum over spins and polarizations, we find 



\M\ 2 = \(Ma + M c )-(Mb + M d )\ 

spin, pol. spin, pol. 



2X2 1 



1 1 

+ ■ 



Tr 



PjlTl^^lPk^Yl^PR 



1 



1 



1 1 

+ 



q\q\ \h-M* Ul -M*J \t 2 - M% u 2 -M* 

2 



Tr 



P/il a ^^i Ph 7V92 Pr + Tr p!iq/2Yl a phqfi l"7 P P R 



+ - 



1 



1 



1 



ql\t 2 -M 2 q U 2 -M* 



Tr 



Pk^Tpkl^l^PR 



(16) 



We evaluate this in terms of scalar products using the 
standard Dirac Algebra, leading to a result too lengthy 
to record here. 

The thermally-averaged rate is given by 



vda=— / ^ X! \M\ 2 dLips 



(17) 



spin, pol. 



where the j arises from averaging over the spins of the 

I 4M 2 

initial \ pair, and v = y 1 is the mean dark matter 

relative velocity, as well as the dark matter single-particle 
velocity in the center of mass frame 3 . 



3 Informative discussions of the meaning of v are given in [36], and, 
including thermal averaging, in [37]. 



The three-body Lorentz Invariant Phase Space is 

dLi p S s = { 2„r d ^ d ^ d3 Q^ p -^- p2 - Q > 



2E\ 2E 2 2Eyy 



(2tt)S 



f!8) 



and P = ki + k 2 . This factorizes into the product of 
two two-body phase space integrals, convolved with an 
integral over the fermion propagator momentum, 



dLips* = f M > d4(^d*p 2 5\P-^ - P2 



'w 2n \2E qi 2E 2 (2tt) 
n d 3 Q 5\ qi -Q- Pl )\ 



2Ei 2E W ' (2tt) 2 J 
J S M2 ^dLips\P\qlpl)dLips\ql,Q\p\)(lV) 



Evaluating the two-body phase space factors in their re- 
spective center of momentum frames, and using P \ — 
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and Mw in the numerator, we find 

1 



g 2 / 4 J , f4 
V(T 512M2 M$tP I x 



In 




4Ag 



1+ln 



- Li 2 



18 



4Af2 



M 2 



(22) 



FIG. 2. The ratio R = va(xx -> e+i/VK")/« a(xx -> e+e") 
for the example model [32], with M% 3> M x . We have used 
« = 10~ 3 c, appropriate for the Galactic halo. 




0.6 0.7 

x w =E w /M x 



FIG. 3. The W spectrum per XX 
ample model, with M x 



evW event for the ex- 

2 ^ A,f2 



300 GeV and M% > M x 



p 2 = 0, we have 



dLips 2 (x 2 ,y 2 ,0) 



y 2 dn 



8lTX 2 47T 

This allow us to write the three-body phase space as 



(20) 



The Spence function (or "dilogarithm" ) is defined as 

U 2 (z) = - /* f In |1 - C| = Er=i W- The ful1 ex P res " 
sion (retaining sub-leading terms in M x in the numer- 
ator) is specified in Appendix D. Clearly, the leading 
terms are neither helicity nor velocity suppressed. 

The effectiveness of the W-strahlung processes in lift- 
ing suppression of the annihilation rate can be seen 
Fig. 2, where we plot the ratio of the W-strahlung 
cross section to that of the lowest order process, Rw = 
V(J {XX e + vW~)/v<r(xx —> e + e~). We see that the 
I'F-strahlung rate rises with DM mass, to quickly dom- 
inate over the lowest order annihilation process. The 
IF-bremsstrahlung rate rises approximately as M*. As 
M 2 increases, eventually phase space allows multi-W/Z 
radiative production, with such a large rate that resum- 
mation techniques become necessary. The onset of multi- 
W/Z dominance has been discussed in [6-8]. 

To obtain the energy spectrum of the W, we compute 
the differential cross section in terms of Ew by making 
the transformation 



dcos(0 q ) -> „, 4 ,^ g 2 dE w . 

(s - q 2 )(q 2 - M w ) 

We find [38], again to leading order M~ 4 , 



(23) 



v da 
dEw 



g 2 f 4 



512E w M w Mfr" 



dLips = 



1 



2 6 (27r)- JM ._ 
{s-qlM-Q 2 ) 



sqf 



dq\ (21) 
d<fr d cos 9pd cos q , 



where <fi is the angle of intersection of the plane defined 
by XX ~^ ee * with that defined by euW, and 0p and q 
are defined in P (CoM) and q rest frames respectively. 

We evaluate the scalar products that arise from Eq.(16) 
in terms of the invariants q 2 , Q 2 — M w , s, t\, and u\, 
and the angles Op, q , and <j>. We then use Eq. (17) to 
evaluate the cross section. As we have neglected dia- 
grams suppressed by M~ 2 relative to those in Fig. 1, we 
present our results to leading order in M~ 4 (i.e., we take 
M 2 3> t±,t2,ui, 112). To leading order in powers of M x 



x ^ 2E W 



E w 



M w (M 2 



w 



6E W + 8E W M X - 2M 2 ) 



AE 



w 



8E W M X 



(2E 



M 2 W ) (2M 2 + M 2 W ) 



x In 



Ew + \J E w 



(24) 



The W spectrum per XX ~^ suW event is given in Fig. 3. 
We use the scaling variable xw = Ew/M x , and plot 



dN/dxw = (j-^) 



da 



dx v 



The kinematic range of 



x w is [j^, (1 + Hp-)], with the lower limit corresponding 
to a IF produced at rest, and the upper limit correspond- 
ing to parallel lepton momenta balancing the opposite W 
momentum. As evident in Fig. 3, the W boson spectrum 
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FIG. 4. The primary lepton spectrum per \X euW for the 
example model, with M x = 300 GeV and M% > M% . 
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FIG. 5. The secondary lepton spectrum (i.e., from W —¥ vgt) 
per W for the example model, with M x = 300 GeV and M% > 
M*. (The branching ratio for W -t vl, 11% per flavor, is not 
included here.) 



has a broad energy distribution, including a significant 
component at high energy Ew ~ M x . 

The energy spectrum of the the primary leptons is cal- 
culated in similar fashion. We present the analytic result 
in Appendix D (along with more detailed expressions for 
v a and vda/dEw)- Here the range of the scaling vari- 
able X£ = Eg/M x is [0, 1 — j&rr]- Both limits arise when 

one lepton has zero energy and the other is produced 
back-to-back with the W. The lepton spectrum is shown 
in Fig. 4. Note that this lepton spectrum is valid for ei- 
ther e + or v from the annihilation \X. ~ * & + vW~ , and 
for either e" or v from the annihilation \X ~~ ^ e~vW + . 
The primary lepton spectrum in Fig. 4 features a sharp 
cut off near Eg = M x , and a dip in the spectrum that is 
due to an absorptive interference effect. 

To obtain the full lepton spectrum, the contributions 
from the subsequent decays of the gauge bosons to lep- 
tons must be included. (The contribution from the lowest 
order 2 — > 2 process XX or vv is negligible. We 



also neglect final state leptons resulting from /x decay and 
from the r decay chain. These leptons are softer than 
those we consider.) For leptons from TU-decay, the range 

of the scaling variable xt is [jjffir, !]• These limits arise 
when all four final state lepton momenta are collinear. 
Particle spectra from the W decay may be calculated 
in a simple but approximate way, as we describe in Ap- 
pendix E leading to Eq. (E8). The resulting secondary 
lepton spectrum is shown in Fig (5). Unsurprisingly, the 
spectrum of secondary leptons is softer than the spec- 
trum of primary leptons. 

When combining the primary lepton and secondary 
lepton spectra, the relative weights are model dependent. 
For example, the primary ^-spectrum is weighted by 
BR{xx -> Wvi)+2BR(xX ZU), while the secondary 
^-spectrum is weighted by BR(\X W+X) x BR(W — ¥ 
vl) + BR(xX -> Z + X) x BR(Z -> U). 

We note that the final charged-lepton spectra will by 
modified by cosmic propagation effects. The injected e^- 
will suffer rapid energy losses from synchrotron and in- 
verse Compton processes on the Universe's background 
magnetic and radiation fields (see, e.g., Ref. [39] for a 
recent analysis) . On the other hand, the injected neutri- 
nos do not interact with the environment, and so their 
spectra remain unmodified. 



C. Unsuppressed Z Emission 

Consider the process producing the vvZ final state. 
The cross sections for the Z-strahlung processes are re- 
lated to those for W-strahlung in a simple way: The am- 
plitudes producing vvZ arise from the same four graphs 
of Fig. (1), where e, W and ~q + are replaced everywhere 
by v and Z and r/o, respectively. The calculation of 
the amplitudes, and their interferences, thus proceeds 
in an identical fashion. After making the replacement 
Mw Mz, the cross section for the annihilation pro- 



cess XX ~^ vvZ differs from that for \X 
only an overall normalization factor, 



-vW~ by 



1 



(2 cos 2 6 W ) 



X V <J e +vW- 



(25) 



~ 0.65 x v o e + vW - 

Consider now the e + e~ Z final state. Again, the am- 
plitudes arise from the same four basic graphs of Fig. (1). 
Since only the left-handed leptons couple to the dark 
matter via the SU(2) doublet rj, only the left handed 
component of e~ participates in the interaction with the 
Z. Therefore, the couplings of the charged leptons to Z 
and W take the same form, up to a normalization con- 
stant. We thus find 



2 (sin 2 0i 



va e+e - z = 



« 3 



cos 2 9 
0.19 x va. 



x v <J e + v w- 



Mi, 



tAf z 



e+vW~ 



Mv 



(26) 
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V. CONCLUSIONS 



In an attempt to explain recent anomalies in cosmic ray 
data in a dark matter framework, various non-standard 
properties have been invoked such as dominant annihila- 
tion to leptons in so-called leptophilic models. When the 
dark matter is Majorana in nature, such annihilations 
invariably are confronted by suppressions of such pro- 
cesses via either p-wave velocity suppression or helicity 
suppression. With the aid of Fierz transformation tech- 
nology, which we have presented in some detail, we have 
elucidated the general circumstances where suppressions 
may be encountered. 

It has been known for some time that photon 
bremsstrahlung may have a dramatic effect on such sup- 
pressions. We have shown that once one considers the 
inclusion of three body final states due to electroweak 
bremsstrahlung, one may also lift these suppressions and 
obtain rates which may be several orders of magnitude 
beyond those without such radiative corrections. In 
fact, barring an unexpected mass-degeneracy, the EW- 
bremsstrahlung lifts the suppression at one order lower 
in a certain small ratio of squared masses than does EM- 
bremsstrahlung, as explained in the text. 

Such radiative processes may be lethal for models 
attempting to produce positrons without overproduc- 
ing antiprotons due to the subsequent hadronization 
of the radiated gauge bosons. Given that electroweak 
bremsstrahlung is the dominant annihilation channel for 
the DM models under consideration, and both W and Z 
decay to hadrons with a branching ratio of approximately 
70%, a large hadronic component is unavoidable. Impor- 
tantly in the context of recent cosmic ray data, there 
will be sizable antiproton production. We also note that 
dark matter searches triggering on anti-deuterons will 
find a sample in the W- and Z-bremsstrahlung processes. 
The Aleph experiment has measured an anti-deuteron 
production rate of 5.9 ± 1.9 x 1CT 6 anti-deuterons per 
hadronic decay of the Z [40]. We expect the rate for 
anti-deuteron production in W-decay to be similar. 

Even for models which do not suffer a suppression of 
the lowest order process, we see that it is impossible 
to have purely leptonic annihilation products, including 
"leptophilic" models in which the dark matter has direct 
couplings only to leptons. In a broader context the re- 
sults presented here show the importance that may be 
played by electroweak bremsstrahlung in future searches 
of indirect dark matter detection. For any DM model for 
which electroweak bremsstrahlung makes an important 
contribution to observable fluxes, there will be large, cor- 
related fluxes of e^, neutrinos, hadrons and gamma rays. 
We will explore the detection of these signals in a future 
article [41]. 
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Appendix A: Fundamentals of Fierzing 

In this paper we have made use of standard Fierz trans- 
formations, helicity-basis Fierz transformations, and gen- 
eralizations of the two. In this Appendix, we derive these 
transformations. The procedure for standard Fierz trans- 
formation can be found in, e.g., [26], while more general 
Fierz transformations are laid down in [25]. The starting 
point is to define a basis {T B } and a dual basis {F^}, 
each spanning 4x4 matrices over the complex number 
field C, such that an orthogonality relation holds. The 
standard Fierz transformation uses the "hermitian" bases 

{T B } = {1,175,7", 757",^"}, and 

{T B } = {1,H75),7^(-757m),7J^}> ( A1 ) 

respectively. Because of their Lorentz and parity trans- 
formation properties, these basis matrices and their duals 
are often labeled as S and S (scalars), P and P (pseu- 
doscalars), V and V (vectors, four for V, four for V), 
A and A (axial vector, four for A, four for A), and T 
and T (antisymmetric tensor, six for T, six for T). As 
usual, spacetime indices are lowered with the Minkowski 
metric, 75 = 7 s = i'y ^ 1 ^ 2 ^ 3 , a^ v = §[7^, 7% (and 
ry5 a ^ — t_ e nvafi Note the change of sign between 
the the basis and dual for the P and A matrices. The 
bases are "hermitian" in that 7 = T B , so that the 

associated Dirac bilinears satisfy [\I>iP B \t'2] t = *2r B *i 
and [^ir B * 2 ] t = * 2 r B *i- Importantly, we have T B = 
(r s ) _1 in the sense of the accompanying orthogonality 
relation 

Tr [T c T B ] =45%, B, C = l,...,16. (A2) 

Note that the factor of ^ in the definition of T = 
(but not in T = a^) provides the normalization required 
by Eq. (A2): 

Tr [T B T B ] (nosum) = Tr Fc r B ] = 4 . (A3) 

c 

The orthogonality relation allows us to expand any 4 x 
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4 complex matrix X in terms of the basis as 
X = X B T B =X B T B , with 
X B = 

1 1 



^Tr[XT B ], and X s = jTr[XT B ] 



i.e., X=^Tr[XT B }T B = -Tr[XT B }T B . (A4) 

One readily finds that the particular matrix element 
(X)ab satisfies 



(X) cd SdbS a 



= -[(X) cd (T B ) dc ](T B ] 



ah ■ 



(A5) 



Since each element (X) cd is arbitrary, Eq. (A5) is equiv- 
alent to a completeness relation 



(i) [i] = \(r B ]\r B ) = \o? B ]\r I 



(A6) 



where we have adopted Takahashi's notation [27] where 
matrix indices are replaced by parentheses (• • • ) and 
brackets [•••], in an obvious way. Thus, any 4x4 matri- 
ces X and Y may be expressed as 



(x)[Y] = (xi) [iY] = i(xr s Y] [v B ) 
= irr[xr B Yr c ] (r c ] [r B ). 



(A7) 



This equation is presented as Eq. (7) in the main text. 
Alternatively, we may express any 4x4 matrices X and 
Y as 

(X) [Y] = (XI) [Yl] = \ (XT B ] [YT B ) (A8) 



1 

43 



Tr [X T 1 



4 



Tr [YT B T U ] (r°] [r 



The RHS's of Eqs. (A7) and (A8) offer two useful op- 
tions for Fierzing matrices. The first option sandwiches 
both LHS matrices into one of the two spinor bilinears, 
and ultimately into a single long trace. The second option 
sandwiches each LHS matrix into a separate spinor bilin- 
ear, and ultimately into separate trace factors. Eq. (A7) 
seems to be more useful than (A8). One use we will make 
of Eq. (A7) will be to express chiral vertices in terms of 
Fierzed standard vertices. But first we reproduce the 
standard Fierz transformation rules by setting X = T D 
and Y = T E in Eq. (A7), to wit: 

(T D ) [T E ] = ±Tr [T D T B T E T c ] (T c ] [T B ) . (A9) 

(An additional minus sign arises if the matrices are sand- 
wiched between anticommuting field operators, rather 
than between Dirac spinors.) Evaluation of the trace 
in Eq. (A9) for the various choices of (B, C) leads to the 
oft-quoted result [26] 



/ (S) [S] \ 
(V) [V] 
(T) [T] 
(A) [A] 
V (P) [P] J 



(l 

4 
6 
4 
VI 



1 

2 


-2 
-1 



1\ 

-4 
6 
-4 

1/ 



/ (S] [§) \ 
(V] [V) 
(T] [f ) 
{A] [A) 

V (P] [P) J 

(A10) 



TABLE II. Fierz-invariant combinations in the standard ba- 
sis. 



Fierz-invariant combination 


eigenvalue 


3(S®S + P(g>P) + T(g>T 


+ 1 


2 (S <g> S - P (8) P) + (V ®V + A® A) 


+ 1 


V ®V - A® A 


-1 


S®S+P®P-T®T 


-1 


2 (S ® S - P ® P) - (V ® V + A® A) 


-1 



More relevant for us, as will be seen, is the ordering 
P, 5, A, V, T, which leads to a Fierz matrix obtained 
from the one above with the swapping of matrix indices 
1^244->3-^54l. The result is 



/ 


1 


1 


-1 


-1 


1\ 


/ (P\ IP) \ 




1 


1 


1 


1 


1 


(S] [S) 




-4 


4 


-2 


2 





(A] [A) 




-4 


4 


2 


-2 





(V] [V) 


V 


6 


6 








-2/ 


V (T] [T) J 



( (P) [P] \ 

(S) [S] 
(A) [A] 
(V) [V] 
V (T) [f] ) 

(AH) 

(The zeroes make it clear that Fierzing induces no cou- 
pling between tensor interactions and vector and axial 
vector interactions.) As an example of how to read this 
matrix, 



(A) [A] = ~(P] [P) + (S] [S) \(A] [A) + \{V] [V), 

(A12) 

or, multiplying by spinors and giving the explicit forms 
of the gamma-matrices, 



(u757 M u) (w(-757p) u ) 
= — (uij5v) (v (— ^75) u) + (uv) (vu) (A13) 

-^("757^) 0(~757m) u ) + \ (u^v) {v-y^u). 



The Fierz matrix M for the standard basis is nonsin- 
gular, and hence has five nonzero eigenvalues Xj. Since 
two swaps of Dirac indices returns the indices to their 
original order, the matrix is idempotent, with M 2 = 1, 
or equivalently, M _1 = M. Accordingly, the five eigen- 
values satisfy A| = 1, so individual eigenvalues must be 
Xj = ±1. Also, the corresponding eigenvectors are in- 
variant under the interchange of two Dirac indices. In 
Table (II) we list the eigenvalues and "Fierz-invariant" 
eigenvectors. 

Helicity projection operators are often present in the- 
ories where the DM couple to the SI/ (2) lepton doublet, 
so it is worth considering Fierz transformations in the 
more convenient chiral basis. 

One derivation of chiral Fierz transformations utilizes 
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the following chiral bases (hatted) [25]: 

{f B } = {Pr, Pl, PsTf, P L r, \°n , and 

{f b} = {Pr, Pl, Pl^, PrJ», \ <V> . ( Al4 ) 

where Pr = 5(1 + 75) and Pl = 5(1 — 75) are the usual 
helicity projectors. The orthogonality property between 
the chiral basis and its dual is 

Tr[t c t B ]=2Sg, B,C = 1,...,W, (A15) 
which implies the normalization 

Tr [V B f B ] (no sum) = ]T Tr t f c T B }=2. (A16) 

c 

Notice that because {75, 7 M } = 0, the dual of Pr7 m is 
Phi ii, and the dual of Pl7^ is PrI^- Notice also that 
the normalization for the chiral bases necessitates factors 
of \ in both T = ■^o >xv and T — \<y^ v , in contrast to the 
tensor elements of the standard bases, given in Eq. (Al). 
In the chiral basis, one is led to a general expansion 

X = 1 Tr [X f B \ ?b = \tt [X f B ] f B , (A17) 

and to a completeness relation 

(!)[!} = 1 -(t B ][t B )= 1 -(t B ][t B ). (A18) 
Thus, any 4x4 matrices X and Y may be expressed as 
(X)[Y] = (Xl)[lY] = i(xf B Y][f B ) 

= ±Tr[Xt c Yt B ] (f B ][t c ). (A19) 

Substituting X = f D and Y = f E into Eq. (A19), one 
gets 

(f D ) [t E \ = l -Tr [f D f c f E t B ] (f B ] [f c ) • (A20) 

Evaluating the trace in Eq. (3) leads to the chiral-basis 
analog of (A10) or (All), presented in Eq. (4) of the 
main text. 



As a check, we note that the matrix M in Eq. (4) is 
idempotent, M 2 = 1, as it must be. The eigenvalues are 
therefore ±1. Eigenvalues and Fierz-invariant eigenvec- 
tors for the chiral basis are given in Table (III) . The final 
two eigenvectors in the Table simply express again the in- 
variance of V ± A interactions under Fierz-transposition 
of Dirac indices. This invariance is also evident in the 
diagonal nature of the bottom two rows of the matrix 
Eq. (4). 

One may instead want the Fierz transformation that 
takes chiral bilinears to standard bilinears. Since mod- 
els are typically formulated in terms of chiral fermions, 
a projection onto standard s-channel bilinears would be 

TABLE III. Fierz-invariant combinations in the chiral basis. 



Fierz-invariant combination 


eigenvalue 


3 (P R <g) Pr + Pl <g> Pl) + f ® f 


+ 1 


2 Pr ®Pl+ Pi?7 M ® PlIu 


+ 1 


2 Pl®Pr + PlY <g> Pa7m 


+ 1 


Pr®Pr + Pl®Pl-T®T 


-1 


2 Pr ®Pl~ Pr J" ® Pl7m 


-1 


2P l ®Pr- PlI 1 " ® PrIh 


-1 


Pfl7 M ® PrIh 


-1 


PljI* ® Pl1„ 


-1 



well- suited for a partial wave analysis. Because differ- 
ent partial waves do not interfere with one another, the 
calculation simplifies in terms of s-channel partial waves. 

Setting X = f D and Y = t E in Eq. (A7), we readily 
get 

(f D ) [t E ] = ±Tr [f D T B t E T C ] (T c ] [T B ) . (A21) 



We (should) get the same result by resolving the RHS 
vector in Eq. (4) into standard basis matrices. The result 
is 
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/ (Pr) [Pr] \ 




/ 1 


1 


1 


1 


1 


1 











o\ 




/ (1] [1) \ 


(Pl) [Pl] 




1 


1 


-1 


-1 


1 


-1 
















(75] [75) 


(p R r) [Pl-tv] 




4 


-4 


4 


-4 






















(75] [1) 






4 


-4 


-4 


4 






















(1] [75) 


(f) [f] 


1 


6 


6 








-2 



















(T] [f ) 


(75T) [f] 


8 








6 


6 





2 
















(75T] [f ) 


(Pr) [Pl] 






















1 


— 1 


1 


-1 




/ in r \ 

(7 M ] [7m) 


(Pl) [Pr] 






















1 


-1 


-1 


1 




(757 M ] [757m) 


{Parf) [PrI A 






















-2 


-2 


-2 


-2 




(7571 [7m) 


\ [PLin [Pl7,] I 


Vo 

















-2 


-2 


2 


2/ 


V m [757m) / 



(A22) 



All relations are invariant under the simultaneous inter- 
changes Pr o Pl and 75 — > —75. The matrix in (A22), 
relating two different bases, is not idcmpotcnt. In fact, 
it is singular. 



Appendix B: Cancellation of Vector and Tensor 
Amplitudes for Majorana Fermions 

Majorana particles are invariants under charge conju- 
gation C. Accordingly, the Majorana field creates and 
annihilates the same particle. This implies that for 
each i-channel diagram, there is an accompanying u- 
channel diagram, obtained by interchanging the momen- 
tum and spin of the two Majorana fermions. The rela- 
tive sign between the t- and u-channcl amplitudes is —1 
in accord with Fermi statistics. For example, consider 
the Fierzed (i.e., s-channel) bilinear for ^-annihilation: 
v(ki, s\)T B u(k2, s 2 )- The associated Fierzed bilinear 
from the (k\ <-» /c 2 )-cxchange graph, with its relative 
minus sign, is — v(k2, S2)T B u(ki, s\). Constraints relat- 
ing the four-component Dirac spinors to their underly- 
ing two-component Majorana spinors must be imposed. 
These constraints, any one of which implies the other 
three, are 



u(p,s) = Cv T (p, s) , u(p,s) 
v(p, s) = Cu T \p, s) , v(p,s) 



-v T (p,s)C-\ 
-« T (p,a)C- 1 .(Bl) 



Here, C is the charge conjugation matrix. These Ma- 
jorana conditions on the spinors allow us to rewrite the 
exchange bilinear as (suppressing spin labels for brevity 
of notation) 

-w(fc 2 )r B u(fc 1 )) = u T (fc 2 )C- 1 r B Ci; T (A :i ) 

= [v{k l ){C- 1 Y B C) T u{k 2 )] T 

= v{k 1 ){n B T B )u{k 2 ). (B2) 

For the final equality, we have used (i) the fact that the 
transpose symbol can be dropped from a number, and (ii) 
the identity {C^TbCV = (vb(T b ) T ) T = r] B T B , where 
rj B = +1 for T = scalar, pseudoscalar, axial vector, and 
r\ B = — 1 for r = vector or tensor. 



In the Four-Fermi or heavy propagator limit, where 
the differing momenta in the t- and u-channel propaga- 
tors can be ignored, one obtains an elegant simplification. 
Subtracting the u-channel amplitude from the t-channel 
amplitude, one arrives at the weighting factor (1 + i] B ), 
which is two for S, P, and A couplings, and zero for V 
and T couplings. Thus, we must drop V and T couplings 
appearing in the Fierzed bilinears of the x _ current. What 
this means for the model under discussion is that after 
Fierzing, only the axial vector coupling of the x-current 
remains, and the factor of 1 + t)a = 2 is multiplied by the 
(7-8)-elemcnt = | in the Fierz matrix of Eq. (4) to give 
a net weight of 1. 



Appendix C: Non-Relativistic and 
Extreme- Relativistic Limits of Fermion Bilinears 

We work in the chiral representation of the Dirac alge- 
bra, and we follow the notation of [28] . Accordingly, 



7o 




7 : 




75 



-1 
1 



The rest-frame four-spinor is 



i(p = 0) 



M 



(Cl) 



(C2) 



where £ is a two-dimensional spinor. The spinor with 
arbitrary momentum is obtained by boosting. One gets 



u(p) 



( Vp^z\ 



(C3) 



where a = (1, (?) and a = (1, —a). 

In a standard fashion, we choose the up and down spin 
eigenstates of cr 3 as the basis for the two-spinors. These 
basis two-spinors are 



(C4) 
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In terms of the chosen basis, we have for the NR u- 
spinors, 



(C5) 



We get the ER limit of the u-spinors from Eq. (C3). After 
a bit of algebra, one finds 




ER 



'2E 



i- 






(C6) 



The arbitrary w-spinor is given by 



v(p) = 



1 VW^V^ 
V -VP ■ a V J 



(C7) 



In the Dirac bilinear the two-spinor rj is independent of 
the two-spinor xi, and so it is given an independent name, 
r\. However, the basis rj^ remains ^ as defined above. 
It is the minus sign in the lower components of v relative 
to the upper components that distinguishes v in eq. (C7) 
from u in eq. (C3) in a fundamental way. After a small 
amount of algebra, one finds the limits 



NR ,^ . 

v± — > v±(p = 0) 




(C8) 



and 




ER 



2E 



(C9) 



Finally, we apply the above to determine the values of 
Dirac bilinears in the NR and ER limits. The u = v) 70 



the u and v spinors. We let T denote any of the hermitian 
basis Dirac- matrices {1, 175, 7^, 75 7^, a^}. Then, the 
NR limit of u(pi) T v(p2) is just 



u( Pl )Tv(p 2 ) ^ M 



(a, 6) r 




(CIO) 



Non-relativistic results for the various choices of basis T's 
and spin combinations are listed in Table I of the text. 

To give a succinct formula for the ER limit of 
u(pi)T v(p2), we take p\ = —p2 = 3, i.e. we work in 
a frame where p\ and p 2 are collinear, and we quantize 
the spin along this collinear axis. The result is 



u(p 1 )Tv(p 2 ) -^4 V4£i^2 



£1 (A+,A_)r 



-A_ 



V2 



(Cll) 

where the matrices A± are just up and down spin pro- 
jectors along the quantization axis 3: 



1 




A_ = 




1 



(C12) 



Extreme-relativistic results for the various choices of ba- 
sis T's and spin combinations arc listed in Table I of the 
text. 



Appendix D: Full Cross Section Results 



We present here the full results of the cross section 



calculations for the process \X 



e^vW^, including 



terms of all orders in M x . In Section IV we presented 
only the leading order terms, which dominate in the large 
M x limit. For M x not too much heavier than Myy, it is 
important to retain sub-leading terms. 

(")x 



and v = 70 conjugate spinors are are easily found from The total cross section for XX C=F v ^ ± is given by 



g 2 / 4 



' 7 M*r 

,32 



7M^ 



-AM^-lGM^Ml-^ 



9 Ml 



Mi 



+ ln 



+ 8M^(M^ + 2M 2 )(Li 2 Ml/Jr \ Ml -Li 2 ^ + %, ) + 0{v\ M~\ m\) 



4M 2 



2M W M X 



-4M 2 



(Dl) 
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The W energy spectrum is 
vd<j c +„ w - g 2 f 4 



dEw 



512E W M^M^7T 3 



2E W JE 2 V - M 2 , (M^r - &E 2 W + 8E W M X - 2M. 



(D2) 



+ (AE§v - %E 3 W M X + (2E 2 W - Mir) (2M X 2 + M 2 



w 



In 



Ew + V E w ~ M w 



Ew - V E w - M w 
while the lepton spectrum (for either the charged lepton or the neutrino) is 

v da c+uW - _ g 2 f 4 \ E e (4(M X - E e )M x - M 2 ,) 



+ 0(v 2 ,M- 2 ,mj) 



dEz 2^M^{M x -E e )^\M 2 v M^{M 2 v +AE e M x )(M x ^E e 



x ( 7 x 2 8 E 7 e M^r + 2 & E & e Ml 



Mlr-^-M 2 ) - '/'E-'M: { ..U 1 



w 



l -fM 2 x M 2 w 



984M- 



2 i E i e M x (Mir - ISMirM 2 - 232M 4 M^ + 2704M^) 



- L ' i -Mir + ^MlrM 2 - SGOM^M 4 - 4672M X 6 M^ + ^^-M* 



1 



1376 



t SE'^M X ^Mlr - 6MlrM 2 + imMlrM* + —M^Mlr - 1600M* 
- 2E e M 2 (M^M 2 - AMlrM 2 + 400M^M 4 + 960M°M^ - 2 10 M*) + 2 & M* (M^ + 2M 2 ) 



- r\2M 2 + M 2 r- 4EA ^ Eef )h, 



M^M X 



(M x -E e )(Mlr + 4E e M x ) 



+ 0(v 2 ,M- 2 ,m 2 ). (D3) 



Appendix E: Approximate Spectrum for Boosted W 
Decay Products 

If any possible polarization of the produced W is ne- 
glected, then a simple calculation results for the spectra 
of the finals state particles from W decay. The lab frame 
spectra of the decay product (of type or "flavor" F) de- 
pends on a one-dimensional convolution of the isotropic 
spectrum in the W rest frame (RF energy £"), ^§r-, with 

the W spectrum in the lab frame, -jp^- We now develop 
this convolution. 

Given the energy distribution dNyy/d-f of produced 
W's (with 7 = Ew/Mw), and the energy distribution 
dNp/dE'p of decay particle F in the W rest frame, 
normalized to the multiplicity of F per W decay (i.e., 
there is a branching ratio W — > F multiplier implicit in 
dNp/dE'p) and assumed to be isotropic, 4 one gets the 
spectrum dNp/dEp of particle F in the lab via: 



dNp(E) 
dE 



1 dcosf?' f , dN w 



(El) 



JdE'^5(E-[ 1 E' + ^p'cos6']), 



4 If the W polarization is not neglected, then the W decay am- 
plitude includes Wigncr functions (9), which introduce a 
linear cos0 or sin# term into Eq. (El). 



with p' = y/E' 2 - m 2 F , /3 7 = x/j 2 - 1- The cos 6' inte- 
gral is easily done, and one gets 



dNp(E) _ 1 
dE 



rf7 dN w f E '+ dE' dN F 



" 2 } Y yff^l d 1 J E , p' dE' ' 

(E2) 

with E'± = "fE ± P"fp. Equivalently, we get 



dNp(E) 
dE 



dE' dN F P+ d 1 dN^ 



w 



p' dE' %_ 



(E3) 

with 7± = (EE' ± pp')/m 2 F and p = x/E 2 -m 2 F . This 
formulation neglects interferences between identical par- 
ticles produced in both the primary and secondary chan- 
nels, if any. 

As given, Eq. (E3) applies to any particle type in the 
W's final state. For example, it could be used to cal- 
culate the antiproton or antineutron spectrum from W 
production and decay, if the fragmentation functions for 
W — > p or n, i.e. f(x§ = 2E§/Mw) were input. 

Here we perform a the convolution for the especially 
simple case of W decay to two massless particles, say v e 
and e. For massless leptons, we have 

( ' X " d -^ i = BR{W^ VC )5{E'- l -M w ) 1 (E4) 



dE' dE' 

with 7+ = (E w /M w ) max = 
(4M% + M$y)/4M X M W , and 7 _ 



(s + M^)/2^M W « 
= (4E 2 + Mlr)/AEM W . 
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The spectrum in the lab is given by Eq. (E3) becomes 
just 



dN v _ dN e _ (BR) 
~dE 



dNe _ (BR) r + 
dE ~ M w J 



dj dNw 
yff^l d 1 



(E5) 



The IF-spectrum shown in Fig. (3) is approximately half 
of an ellipse, suggesting the fit 



In 



dN 
d:r\\ 



■)- 



In 0.07 



In2.0-ln0.07 



xw — 0.65 
0.50 (1.01 - 0.29) 



= 1 . 



(E6) 



valid for 0.29 < x w < 1.01. Solving for dN/dx w then 
gives 



dN ( 2 \ V 1 - 7 - 7 ^"-" - 65 ) 2 

= 0.07 1 



dxw 



0.07 



(E7) 



M w dl 



M w dN 
dx\ V 



Substituting into Eq. (E5) 7 

and j^- given in Eq. (E7), we obtain the desired one- 
dimensional integral for the secondary lepton spectrum, 
per W: 



dN„(x e ) 
dxg 



= 0.07 (BR) / 

J X- 



dxw 



x w 



2.0 
O07 



{x w -0.&'o)' 2 



(E8) 



The integration limits are x + 



EL 

4M3 ' 



(x£ + J ■ The range of xg for the leptons from 

W decay is [jjp-, 1]. The relevant branching ratios [42] 

are BR(W -> ve) = 11%, BR(Z — > vv) = 6.7%, and 
BR(Z — » t + l~) = 3.4%, each per single flavor mode, 
e, fi, or t. We show the resulting lepton spectrum, with- 
out the BR factor, in Fig. (5). 
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